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Abstract. We compare different approximation schemes for investigating ferromagnetism in the periodic 
Anderson model. The use of several approximations allows for a detailed analysis of the implications of the 
respective methods, and also of the mechanisms driving the ferromagnetic transition. For the Kondo limit, 
our results confirm a previously proposed mechanism leading to ferromagnetic order, namely an RKKY 
exchange mediated via the formation of Kondo screening clouds in the conduction band. The contrary case 
is found in the intermediate-valence regime. Here, the bandshift correction ensuring a correct high-energy 
expansion of the self-energy is essential. Inclusion of damping effects reduces stability of the ferromagnetic 
phase. 

PACS. 71.10.Fd Lattice fermion models (Hubbard model, etc.) - 71.28.-)~d Narrow-band systems; 
intermediate-valence solids - 75. 30. Mb Valence fluctuation, Kondo lattice, and heavy-fermion 



1 Introduction 

The periodic Anderson model (PAM) provides probably 
the simplest starting point to investigate the interplay be- 
tween the hybridization of a periodic array of localized 
(/-) electron levels and a band of itinerant electrons, the 
strong correlation of the localized electrons and quantum 

■ mechanics (Pauli principle) [Q . 

Within this model, we will investigate the many-body 
problem of ferromagnetic ordering. Ferromagnetism in the 
PAM has previously been examined, but most efforts were 

■ aimed solely at the so-called Kondo regime of the model. 
' This is defined by an integer number of /-electrons per 
, lattice site. These localized electrons can form an array 

of well-defined local moments. Via an RKKY interaction, 
' these can order ferro- or antiferromag netically|,|,|||. 

■ In previous papers, we have shown that in the 
intermediate- valence regime of the PAM, defined by a 
non-integer filling of the localized levels accompanied by 
a larger effective hybridization, ferromagnetic order is 
also possible[^^||,^. This, however, raises some ques- 
tions concerning the driving mechanism for the ferromag- 
netic ordering. First of all, in the intermediate-valence 
regime, quantum fluctuations effectively delocalize the 
/-electrons, thus the formation of stable moments needs 
further explanation. A second question is due to the fact, 
that some of the methods used in the above-mentioned 
papers, exclude by construction such effects as RKKY 
as we will discuss below. So why can these methods 
give meaningful results when neglecting something con- 



sidered as essential as RKKY exchange? What should 
be the driving force towards ferromagnetic ordering in 
the intermediate valence regime? In this paper, we try 
to clarify the apparent incoherence of the investigations 
of ferromagnetism in the PAM in the Kondo and the 
intermediate-valence regime. 



After introducing the model in section 2.1, we will dis- 
cuss sev eral well-known propert ies of the system in sec- 
tion 2^ and finally, in section 2_^, introduce a number 
of approximation schemes. By comparing to the known 
properties discussed before, the advantages and disadvan- 
tages of these methods will become clear. In section ||, 
we present and compare results obtained with the differ- 
ent methods. Knowing the strengths and shortcomings of 
them will help us to understand the origin of ferromag- 
netic order. Altogether, this will lead us to the conclusion 
that there are indeed two distinct mechanisms at work: in 
the Kondo regime, an RKKY exchange has to be seen as 
the cause of the ferromagnetic order. In the intermediate- 
valence regime, however, the situation resembles more that 
of a band-ferromagnet as described by the single-band 
Hubbard model. In section ^ we will conclude with a sum- 
mary of our findings. 



2 Theory 

2.1 The Periodic Anderson Model 

We investigate the standard form of the periodic Ander- 
son model (PAM) where a non-degenerate localized /-level 
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hybridizes with a non-degenerate conduction band ("s- The spin-dependent average occupation number 7^^*''^'' can 
band") via an on-site hybridization: now easily be determined: 



k.cr i,a 

i,(7 i.{T 

Here, Sko- [fia) and s\^^ {.fja) are the creation and anni- 
hilation operators for a conduction electron with Bloch 
vector k and spin a (a localized electron on site i and spin 

a) and n-^J = fl^f^a- Ska = :^ Ek e'''^'Si<T and e(k) is the 
dispersion of the conduction band and Cf is the position of 
the localized level. The hybridization strength V is taken 
to be k-independent, and finally, U is the on-site Coulomb 
interaction strength between two /-electrons. Throughout 
this paper, the conduction band will be described by a 
Bloch (free) density of states, po{E) — J2k ^(-^ ~ ^C*))' 
of semi-elliptic shape. Its width W = 1 sets the energy 
scale, and its center of gravity the energy-zero: Tu = 

wi:k<k) = o. 

The problem of determining the relevant (Zubarev) 
Green functions JlO|, pi] , 



G 



N ^ 



I'ja 



{E) 



(2) 
(3) 



can be reduced to the a priori unknown self-energy via 
the formal solultion of the respective equations of motion: 



g\1}{e) = 



E - {ei - fi) - EkaiE) 



{E - (ef - ^i) - Ek4E))iE - (e(k) - - 

(4a) 

1 



E-ia- fi) - E-(e(U)-^) 

Here, the self-energy is defined by 



(4b) 



SME)Gi^]iE) = C/1 5]((/^_./q-./p + k-q.; fL) 



(5) 



Throughout this paper, we will apply the local approxi- 
mation, i. e. assume a k-independent self-energy. Although 
becoming exact only in the limit of infinite spatial dimen- 
sions, it was shown that this approximation gives satisfac- 
tory results already for three dimensions [ L2i[l^ . 

From the Green functions ( |4a| ) and (4^), the /- and 
s-quasiparticle densities of states (/- and s-DOS) can be 
calculated: 



is) 



iO^ 



k 



(6a) 
(6b) 



/ + 00 
dEf_{E)pi^\E) (7a) 
-OO 

/ + 00 
dEUE)pif\E) (7b) 
-oo 

Here, f-{E) denotes the Fermi function and 3x the imag- 
inary part of X. 

Before introducing our approximative methods, let us 
discuss some general properties of the PAM. 



2.2 General Properties 

2.2.1 The Hybridization-Free Case 

For vanishing hybridization strength V ( "atomic limit" ) , 
the problem reduces to that of the zero-bandwidth Hub- 
bard model llj, [ill . One obtains for the /-Greenfunction: 



g,(/,at.) ^ 



E^{e,-p)- U{1 ~ n^ll) 



{E-{ei-p)){E-{ei-p)-U) 
with the respective self-energy: 

Un^^liE - (et - fi)) 



sI^^-\e) 



E - (cf - /i) -U{\-n 



(8) 



(9) 



The excitation spectrum consists of two peaks located at 
Cf and et -I- U , which are called charge excitations. The 
conduction band DOS remains unchanged. 



2.2.2 The Non-Interacting Limit 

The second trivial limit, the interaction- free limit {U — 0), 
yields the following /-Green function: 



G'^r'\E)= 



1 



E-{ef- p) 



B-(e(k)-M) 



(10) 



The DOS now consists of two features, one correspond- 
ing to the /-level, which becomes broadened due to the 
hybridization. The other feature is the renormalized con- 
duction band. Due to the hybridization there is also an 
admixture of /-spectral weight into the conduction band 
region and vice versa. The amount of this admixture of 
/- and s-states can be understood as a rough estimate of 
the effective hybridization. It is generally stronger for Cf 
close to, or inside the conduction band as for ef well be- 
low the band. A further effect is clearly visible when Cf 
lies within the conduction band. Level-repulsion between 
/-level and the conduction band induces a gap located ap- 
proximately at ef . For ef below the conduction band, the 
level repulsion appears only in form of a small shift of the 
lower edge of the conduction band and the center of the 
/-level is situated slightly below ef. 
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Fig. 1. Density of states for V = 0.2, et = -0.75 and -0.25 in 
the interaction- free hmit {U = 0). Thick hnes: /-contribution, 
thin hnes conduction band contribution. 



These two situations are clearly distinct, one expects 
different findings in each. The situation of ef well below the 
conduction band leads for finite (large) U to the Kondo 
regime of the PAM. Here the /-level is almost integer- 
filled. The opposite case, with Cf located within the band 
and a hybridization strong enough to lead to non-integer 
•n!'^^ is called intermediate- valence regime. Both situations 
are plotted in figure |^. 



2.2.3 High-Energy Expansion 

Next, we introduce a useful high-energy expansion for the 
/-electron Green function and the self-energy: 

oo ,,,(n) 



•' ri=0 



oo ^(ri) 
n=Q 



(11) 



Ska{E) — —^^G^JiE) is the spectral density. Its mo- 
ments, which are the coefficients of the high-energy ex- 
pansion of G^^ [E) , are defined by 



12) 



M^i^J = / dEE'^SuAE); n = 0, 1, 2, . . . 
and can be calculated independently via 

MiZ^ = {[[.. .[h,,H]^,...,H]^jl] + ) (13) 



n-fold commutator 



where [...,...]_ denotes the commutator and]. ..,...] + 
the anti-commutator. By inserting ( |ll|) into (|4b|) , one can 
determine the coefficients of the self-energy expansion. 
For the approximative approaches presented below, we 



coefficients = ^ Ek C^'^' 
A/(o) =1 



M(i) =ef + Uri^ll 



(14a) 
(14b) 

M^^) ^ej + 2eiUn^ll + U'^n^ll + V'^ (14c) 

M^3) + Ze^fUn^Il + t/^ej (2nL^] + nL^]') + U^r^l + 

+ l^'(2ef + 2Un[!l + Tu) + U^n^Il{l - n^Il)B^, 

(14d) 

C(°) = Un^H (15a) 
CW =C/2nL^](l-nL^]) (15b) 
C7(2) = U^n^Il(l - n^ll){B^, + C/(l - n^H)) (15c) 



and 



The abbreviation B^^r in (^_4|) and (|15|) stands for a higher 
correlation function called bandshift: 

n^Ilil - n^ll){B^, ef) = V {fl,s,^,{2n\fj - 1)) 

(16) 

In spite of the fact that it is a " higher" correlation function 
it can rigorously be expressed by the Green function ( ^ ) 
and the self-energy (p[)p5[|: 

n^!l{l~n^!l){B^,-e,) = 
-l^l " dEf_{E) i^^^aiE) - l) (17) 

x({E-{e,-^,)-S,{E))G\il{E)-l) 

Surprisingly the hybridization V does not explicitely ap- 
pear in the C^"'' . The contributions via the moments ( p^ ) 
are exactly cancelled by those from the term ^_^^^^y^^_^\^ 
in (4b). From this, one conclusion can already be drawn: 



only need the local moments Mi-"-* = M'^J and 



An) 



The hybridization V enters the calculation only via equa- 
tion (4b) in combination with the free conduction band 
dispersion e(k). Although there are clearly correlation- 
induced effects in the conduction band (cf. equation (Q)), 
these do not feed back into the determination of the high- 
energy features of the self-energy. Any RKKY-like indirect 
exchange between /-sites driven by correlations (cf. refer- 
ence 1^), leaves no footprints in the high-energy behaviour 
of the self-energy. Although this reasoning works only in 
the local approximation, it does not imply that the local 
approximation itself suppresses any RKK Y exchange (cf. 
discussion and references in section 2.3.4). 

The prominent high-energy features of the PAM are 
the charge excitations known from the zero-hybridization 
limit. One can easily check that in this limit the mo- 
ments ( |l^ ) and self-energy coefficients (|l^) are fulfilled 
since i?_o- — > et. So obviously, the bandshift, and there- 
fore the n = 3-moment take care of a correction of the 
positions and weights of the charge excitations in the case 
of finite interaction and hybridization. We will show below 
that this correction can be decisive for a proper descrip- 
tion of ferromagnetism in the PAM. 
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Fig. 2. /-Density of states as obtained within the MPT (see 
text) for V = 0.2, (7 = 2, n'-'^' = 1.0 and different n*"' as 
indicated, e/ was determined such that rS^^ — 1 holds. The 
chemical potential is positioned at the arrows. The inset shows 
the full energy range, the main picture just the region around 



2.2.4 Low-Energy Properties 

The PAM is the lattice-periodic extension of the single- 
impurity Anderson model (SIAM). Since the latter is 
famous for its special low-energy features ("Kondo- 
physics" ) , one expects similar findings also for the 
PAM. The most prominent finding for the SIAM is the 
Kondo-screening: At low temperatures, the magnetic 
moment of the impurity-site is screened by conduction 
electrons. The remaining conduction electrons form a 
Fermi liquid. All physical quantities can be scaled by 
a single energy: the Kondo temperature Tk- In the ex- 
citation spectra, the most significant signature of the 
screening is the occurence of a sharp resonance at the 
chemical potential, the Kondo Resonance. 

The low-energy properties of the PAM have been the 
subject of intense research||l|,^|l^,0,|l8|,^. For the sym- 
metric PAM, defined by n") ^ 1.0 and Cf = 
a Kondo resonance appears centered at /x. But contrary 
to the SIAM, it is split by a gap, the coherence gap 
and the system is insula ting . Picturing the Kondo res- 
onance as virtual /-level po], the gap is simply due to 
level-repulsion between the "flat band" of the virtual 
/-levels at every lattice-site and the conduction band. 
Although it originates from the same mechanism as the 
gap discussed in section 2.2.2, it is clearly distinguishable: 
The coherence gap is, together with the Kondo resonance, 
pinned at /i whereas the hybridization gap discussed in 
the interaction- free limit would show up at et. In figure H, 
we present the densities of states as obtained with the 
below-discussed modified perturbation theory (MPT). 
Although these DOS are calculated within an approxima- 
tive method, the qualitatively same picture emerges from 
the numerically exact Quantum Monte Carlo (QMC)[§, 
pi| and numerical renormalization group calculations 
(NRG) Hi and can thus be believed to be qualitatively 
correct. 

The situation of the symmetric PAM is plotted at 
the bottom of figure ^. Moving away from the symmetric 



PAM, e.g. by reducing the number of conduction electrons, 
the system becomes metallic again. The Kondo resonance 
gets asymmetric relative to /i, the coherence gap is shifted 
away from the chemical potential. Further away from the 
symmetric parameters, the coherence gap closes, which is 
probably due to quasiparticle damping, since '^Ea{E) ~ 

A still open question is that of a unique energy scale 
similar to Tk for the SIAM. The main problem un- 
der discussion is the exhaustion problem introduced by 
Nozieres||2^. For the case of a periodic array of localized 
/-levels, the system cannot make available enough con- 
duction electrons to screen all /-moments. This situation 
is still subject to many investigations H, |l|,|l8[|l|,||. 



2.3 Approximative Solution of the PAM 

In the following sections, we will introduce a series of 
approximation methods to determine the self-energy (|^) 
of the periodic Anderson model. This series represents 
a subsequent improvement of the theory along the lines 
discussed in the last section: The Hubbard-I approxima- 
tion recovers the non-interacting and the hybridization- 
free limit but fails to reproduce the high-energy behaviour 
of the self-energy. This is corrected by the spectral den- 
sity approximation. This method, however, still suffers 
from the complete neglection of quasiparticle damping. 
A straightforward method to incorporate this is the mod- 
ified alloy analogy. Finally, the low-energy behaviour can 
be qualitatively reproduced by the modified perturbation 
theory. 



2.3.1 The Hubbard-I Approximation 

The first approach presented here is a crude interpolation 
between th e two exactly solvable limitin g case s V = Q 
(cf. section |2.2.lD and ?7 = (cf. section The re- 

sult equals that of Hubbard's first work on the Hubbard- 
model||lJ|, the Hubbard-I approximation. 



In the interaction- free limit, one can express the /- 
Green function ( p^ in terms of the corresponding [7 = 0- 
V = { "atomic" ) solution: 



G 



(f,v=o) 

kcr 



{g9 



(/,U=0,at.) 



(18) 



£;-(e(k)-p) 



The Hubbard-I approximation is now obtained by assum- 
ing the functional dependence of ( p^ also for the finite- 
U case. With the atomic-limit Green function for the full 
PAM (^) , this essentially corresponds to inserting the self- 
energy (ph into equation (4b), therefore 



Un^lljE ~ jet ~ fi)) 



E - (ef - - Uil - n"!^) 



(19) 



Although this method is by construction exact in two 
limiting cases, namely the interaction-free (U ~ 0) and 
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the hybridization-free (V — 0) case, a number of short- 
comings foUows directly from inspecting equation ( |l9| ) and 
its derivation. First of all, the self-energy fulfills the 
high-energy expansion only up to the n = 1 self-energy co- 
efficient. The bandshift correction B-o- is neglected. This 
yields a weak spin-dependence through the expectation 
values n^Jl — {n^}^) only. Second, the self-energy 
is real, any quasiparticle damping effects as indicated by 
a finite imaginary part of S„{E) are ignored. The third 
drawback is the suppression of any feedback mechanism 
from the conduction band: The latter enters the calcula- 
tion only via equation This implies that although 
there are correlation-induced changes in the conduction 
band (cf. equation (|4a|)), these do not feedback into the /- 
self-energy (p^. An indirect magnetic exchange between 
the /-electrons via polarization of the conduction band 
("RKKY") cannot appear in this approximation method. 
Finally, none of the expecte d low -energy features of the 
PAM, as discussed in section 2.2.4 can be found using this 
method. 



2.3.2 The Spectral Density Approximation 

The spectral density approximation (SDA)||2j,|ll|,^ is 
the result of a direct improvement of the above-discussed 
Hubbard-f approximation with respect to the high-energy 
expansion (pT|). 

Starting with an ansatz for the self-energy using the 
same functional structure as (fol). 



E - Q!2ct 

E - aza 



(20) 



one can fit the coefficients ap„ in such a wa; 



that the 

high energy expansion of the self-energy ([flh with the 
coefhcients dl^) is fulfilled. One readily arrives at 



^(SDA)(^) 



Un'^Il{E-B-,-{ei-^i)) 
E-B^,- (ef - /i) - U{1 - n'-Il) 



(21) 



The SDA self-energy differs from the Hubbard-I solution 
by the ba ndshi ft i?_cr ([T^). It is introduced by the n = 3- 
moment (14d) and responsible to reproduce the correct 
high-energy behaviour of the /-Green function. It leads 
to a (possibly spin-dependent) shift of the positions of 
the /-peaks in the density of states. So without loosing 
any of the advantages of the Hubbard-I approximation, 
as the correct reproduction of the U — and the V — 
limits and its numerical simplicity, one major drawback 
of this method can be removed. However, the other points 
of criticism, as the missing quasiparticle damping and the 
incorrect low-energy properties remain. 

The name "spectral density approximation" stems 
from its application to the Hubbard model, where this 
approach is derived by a physically motivated two-pole 
ansatz for the spectral density|^,^. Results for the SDA 
in the context of the PAM were previously published 
In these papers, the SDA was applied to a set of effective 
Hubbard models, onto which the PAM could be mapped. 



However, the above-described procedure leads to _exactly 
the same results as those published in references 



2.3.3 The Modified Alloy Analogy 

Now we want to present a method that resolves one ma- 
jor drawback of the already introduced methods: Both the 
Hubbard-I and the SDA self-energies are real. Quasipar- 
ticle damping as represented by a finite imaginary part of 
the self-energy is therefore completely neglected. A well- 
known method to incorporate damping effects is the alloy- 
analogy approach. By using physical intuition or by other 
justified means, the original problem is mapped onto an 
fictitious alloy, which subsequently can be solved using 
standard methods as e. g. the coherent potential approx- 
imation (CPA)|6|. 

The CPA represents the best "single-site" method for 
solving an alloy-problem Single-site approximation in 
this context is equivalent to the already introduced local 
approximation or k-independence of the self-energy. The 
alloy is defined by the energy levels of its components, Ep^ 
and their respective concentrations Xpa , where the index p 
numbers the components. The corresponding self-energy 
can be determined by solving the CPA-equation: 



= 



E.^ 



pa 



Cf 



1 - Gi{}iE){Ep, - S^iE) ~ e,) 



(22) 



In the conventional alloy analogy for the PAM p7|, Pq, p9[ , 
the artificial alloy is determined by the poles of the atomic 
limit /-Green function (0) and their respective weights: 



2a 



ef H 



JAA) 

''1(7 

'^2a 



1 — n 

if) 



if) 



(23) 



n 



This choice, however, is in no way predetermined. In refer- 
ence M , another alloy analogy, the modified alloy analogy 
(MAA) was proposed for the PAM in analogy to the MAA 
for the Hubbard model|3^: By inserting equation ([ll| ) 
into ( |2^ ) and comparing the coefficients in an opti- 
mum two-component alloy analogy with respect to the 
high-energy behaviour can be found: 



El,2a + C/ + ef± 



(24) 



(B_, + U- ef)2 + AUn^Hiei - B_^)] 

E2a - Cf - Un^Il ^ 
Xlo- = = 1 - X2a 

^2a — -C'lcr 

It should be noted that Ei^2a coincide with the poles of the 
SDA Green function (cf. equations ( pi| ) and (4b)) when 
the term — — nrrr; — -T in the denominator of the Green func- 



E-(e{k)-f,) 

tion is neglected. 

By construction, the MAA fulfills the high-energy ex- 
pansion ( pT| ) up to the n = 3 moment or equivalently 
the n = 2 self-energy coefficient (15c). Also, damping ef- 
fects are considered. However, the MAA still suffers from 
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the co mplete neglection of the low-energy properties (cf. 



section 2.2.4 ) and the self-energy is determined in such a 
way that correlation effects within the conduction band 
do have no influence ("feedback") onto IJa{E). A more 
detailed discussion of the MAA applied to the PAM can 
be found in reference B. 



2.3.4 The Modified Perturbation Theory 

Let us finally introduce the modified perturbation theory. 
This approach is based on the dynamical mean-field the- 
ory (DMFT)(3ll|3|l. In the case of a local self-energy, 
which becomes exact in the limit of infinite spatial 
dimensions [ ^3| , |3^ , the PAM can be mapped onto a single- 
impurity Anderson model (SIAM) with the Hamiltonian, 

H = ^(e(k) - Ai)4,sk. + - M)4rf. + (25) 

k,cr a 

a (7 

The notation is as for the PAM (cf. equation SJ'' are 
the annihilation (creation) operators for electrons at the 
impurity, its energy level is Cd- In the case of the DMFT, 
the bath of conduction electrons, usually defined by e(k) 
and Vkd need not be specified in detail. For all practi- 
cal calculations, it is sufficient to know the hybridization 
function A{E), in the pure SIAM defined by A{E) = 

X^k g-e(k) • However, for the mapping of the DMFT to 
be successfull, the hybridization function has to be deter- 
mined according to the self- consistency condition ^^ 



(26) 



This implies that in case of symmetry breaking, the hy- 
bridization function becomes spin-dependent. Now one 
can make use of the fact that the impurity self-energy of 
the SIAM defined by equation (26) is equivalent to the 
self-energy of the PAM. The advantage of the mapping 
is that the SIAM is one of the simplest know many-body 
models, several exact statements as well as well-tested 
approximative solutions are known. In the following, 
we apply the modified perturbation theory, which was 
explained and discussed in detail elsewhere |^ , 36 , 19|, p| . 
So we will restrict ourselves to a short summary of this 
approach. 

Starti ng poi nt is the following ansatz for the self- 



energy y 



(SOC) 



1 - 



(SOC) 



{E) 



(27) 



aa and (3a- are introduced as parameters to be determined 
later. E^'^'~^\e) is the second-order contribution topcr- 
turbation theory around the Hartree-Fock solution pO[,pl|, 
. Equation ( P7|) can be understood as the simplest pos- 
sible ansatz which can, on the one hand, reproduce the 



perturbational result in the limit [/ — > 0, and, on the other 
hand, recovers the atomic limit for appropriately chosen 
aa and 

Using the perturbation theory around the Hartree- 
Fock solution introduces an ambiguity into the calcula- 
tion. Within the self-consistent Hartree-Fock calculation, 
one can either choose the chemical potential to be equiva- 
lent to the chemical potential of the full MPT calculation, 
or take it as parameter fi to be fitted to another phys- 
ically motivated constraint. In reference |^ and other 
papers |^,|44j, the Luttinger theorem^^ or equivalently 
the Friedel sum 

rule||,|3 

was used to determine ft. Since 
these theorems are applicable only for T = 0, this limits 
the calculations to zero temperature. In order to access fi- 
nite temperatures, we used the condition of identical elec- 
tron densities for the Hartree-Fock and the full calcula- 
tion {ua ' = Ua ). In our view, it is more reasonable to 
perform the Hartree-Fock calculation for the same elec- 
tron density rather than identical chemical potential of 
the full MPT calculation since the electron density is a 
critical parameter concering correlation effects. A more 
detailed analysis of the different possibilities to determine 
fi is found in reference | ]35| . Finally, the parameters and 
j3a have to be determined. Instead of using the "atomic" 
limit of y = as was done e. g. in references [^^, 
we make use of the moments of the spectral density. 
Analogously to equations (pT|), ( p^ and (p^), these can 
be evaluated for the SIAM. To fit the two parameters 
of ansatz (p7|), the first three self-energy coefficients are 

needed, since ci°^ is reproduced for any choice of and 
P.- 

As for the PAM, a bandshift correlation function sim- 
ilar to ( |l6| ) is introduced via C^^ and the procedure leads 
to the correct high-energy behaviour of the Green func- 
tion for the SIAM and via the DMFT-mapping also for 
the PAM. So while recovering the main advantage of the 
SDA and MAA, namely the correct reproduction of the 
high-energy expansion ( [ll| ) up the n — 3-moment, the 
MPT yields a major improvement concerning the low- 
energy properties of the PAM. Although already for the 
SIAM the low-energy scale ("Kondo temperature" Tk) 
connected with these properties, cannot be quantitatively 
reproduced, other quantities caii. at least in a qualitatively 
satisfactory way, be recovered |3^,^. In particular, the 
densities of states both above and below the Kondo tem- 
perature, but also the general features of the susceptibility 
x{T) seem to be trustworthy. Another test of the low- 
energy properties is given by the Friedel sum rule, which 
links the self-energy at the Fermi energy with the electron 
density. Within the MPT, it is fulfilled in a large param- 
eter space p^. It is also worth mentioning that via the 
DMFT self-consistency (|2^), a feedback from correlation- 
induced features in the conduction band onto the /-self- 
energy is possible. In the limit of infinite spatial dimension 
(d — > oo), where the DMFT becomes exact, the RKKY ex- 
change between any two lattice sites will vanish. However, 
as discussed in reference [|| , the net exchange of one lat- 
tice site with a shell of neighburs remains finite since the 
number of sites in the respective shell diverges as c? — » oo. 
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Fig. 3. s-DOS (/-DOS) as dotted (solid) lines for V = 0.1, (7 = 1 and n''°') = 2.0. The respective value of e/ is given in 
each graph. The position of the chemical potential is denoted by the arrows. The left column was calculated with the SDA, the 
middle within MAA and the right column within MPT. 



This exchange is incorporated in the (spin-dependent) hy- 
bridization function (|2q). 



3 Results and Discussion 

In the following we will present the results obtai ned 
with the different approximation schemes of section ^.3| 
and hope to shed some light on the mechanism that 
leads to ferromagnetism in the PAM in the Kondo- and 
intermediate-valence regime. First however, we will look 
at the paramagnetic quasiparticle densities of states 
(DOS) as defined by equations ( pa| ) and (6b). 

In figure ^, both the /- and s-DOS are plotted for a 
relatively small interaction strength C/ = 1.0 and V = 0.2, 
^(tot) _ 2.0 at zero temperature and for various Cf as 
indicated. The second picture from the bottom represents 
the above-introduced symmetric case with ef = —0.5 = 



— ■J. The position of the chemical potential /i is indicated 
by the arrows. The left column was obtained using the 
SDA, the middle MAA and the right column by using the 
MPT. 

The SDA DOS differs from the interaction- free case (cf. 
section 2.2. 2| ) plotted in figure |l| by the appearance of a 
second charge excitation approximatly a,t Cf ^-U .11 either 
Cf ov -\- U falls within the band region, a hybridization 
gap as discussed in section 2.2.2 is clearly visible. The 



Hubbard-I results are not shown in figure |3[ These look 
very similar to the SDA DOS, only a small shift of the 
charge excitations can be noticed. In the symmetic case 
they are identical, the bandshift vanishes. 

The MAA DOS shows some modifications when com- 
pared to the SDA: The quasiparticle damping softens the 
charge excitations and the hybridization gap is for ef < 
almost closed. For ef > 0.25 the DOS strongly resembles 
the SDA results. This can be understood since the num- 



8 



D. Meyer, W. Nolting: Ferromagnetism in the Periodic Anderson Model 




0.2 





■ — ■ H-I 




SDA 




. MAA 




C\ MPT 










Kl-^^.—U ' 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



0.5 



1 

(tot) 



1.5 



Fig. 4. s- and 

et = -0.35 (s: 
I, SDA, MAA 
scales of the y- 



/- polarization for U 
lower, /: upper panel) 
and MPT are shown. 



= 4:,V = 0.1, T = and 
The results of Hubbard- 
Please note the different 




- 1 ' 

PM 

y 


/iFM 

/ 


■FM J 






PM 


- 1 , 


1 , - 




Fig. 5. r = phase diagram for the PAM in the intermediate- 
valence regime with U = 4 and V = 0.2 as obtained within 
SDA, MAA and MPT calculations. For the MAA, the phase 
diagram for 1/ = 0.1 is additionally plotted as thin lines. The 
dotted lines separate the regions of parallel (below the dot- 
ted line) and anti-parallel (above) s-f coupling. Note: for the 
MAA with V = 0.2, the low-density ferromagnetic region ex- 
hibits only parallel, the high-density region only antiparallel 
s- /-coupling. 



ber of /-electrons, n*^-^^ is very small. Scattering processes 
become rare and the quasiparticle damping, which dif- 
ferentiates between SDA and MAA, neghgible. As in the 
case of the SDA and Hubbard-I approximation, the sim- 
pler theory disrespecting the high-energy expansion and 
neglecting the bandshift correction, in this case the con- 
ventional alloy analogy (^3|) , yields very similar and in the 
symmetric case identical results. We will see below that 
the bandshift correction becomes much more important in 
the ferromagnetic phase. 

Finally, the DOS obtained by MPT shows remark- 
able differences, especially close to the chemical poten- 
tial. T hese r epresent the "Kondo physics" discussed in 
section 2.2.4. Again, we note that for Cf > 0.25 the DOS 



resemble the SDA and MAA results. This is obviously for 
the same reasons as discussed above, namely the small 
number of /-electrons. 

In figure ^, the /- and s-magnetization is plotted as 
function of the total electron density n*^*"*^ for U — A, 
a = -0.35, V = 0.1 and T = 0. Within the Hubbard-I 
approximation, a small region of ferromagnetism is found 
around n^*"*^ « 1. The conduction band magnetization 
m*^*^ (thin lines in figure ^) is always positive. We call this 
situation parallel s-f coupling. 

In the SDA the region of ferromagnetism is strongly en- 
larged. This is a clear indication of the importance of the 
bandshift correction with respect to ferromagnetism. Since 
the bandshift (^6|) can be spin-dependent, it enhances the 
possibility of ferromagnetic ordering. We also note an in- 
teresting behaviour of the conduction band magnetization: 
as function of electron density, it changes sign. For low 



^(tot)^ the s-f coupling is parallel, for higher values an- 
tiparallel. This behaviour, which is also found within the 
MAA and MPT, can be traced back to the appea rance of 
the hybridization gap discussed in section 2.2.2, A more 
detailed investigation of this can be found in references 

i- 

The MAA result can be used to estimate the influence 
of quasiparticle damping. One observes a reduction of the 
ferromagnetic region compared to the SDA result. So, sim- 
ilar to the Hubbard model[^0[^, quasiparticle damping 
is unfavourable for ferromagnetism. However, apart from 
the lower critical the MAA and SDA curves are 

very similar. It should be noted that with the conven- 
tional alloy analogy ( p3| ) , no ferromagnetic solution can 
be found p^]. This confirms again the importance of the 
bandshift correction ( |l6| ) with respect to ferromagnetism. 

The same holds true for the MPT result. Again, the 
region of ferromagnetism is reduced as compared to the 
MAA, but the change is rather small. Qualitative features 
as the change of sign of the conduction band magnetiza- 
tion, and the generally larger conduction band polariza- 
tion in case of parallel coupling remain the same for SDA, 
MAA and MPT. It should further be pointed out that the 
lower critical n(*°'' is, when varying Cf , in fact determined 
by n'-^\ The number of conduction band electrons, n^"' 
plays no significant role|9|. 

From all these observations we conclude that the band- 
shift correction has a strong influence on ferromagnetism. 
So does the inclusion of quasiparticle damping. However, 
at least in the examined situation of the PAM in the 
intermediate-valence regime, there seems to be no ma- 
jor difference between the MAA without, and the MPT 
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Fig. 6. The Curie temperature as function of total electron 
density for (7 = 4, V = 0.1 and et = -0.35. 



including the special low-energy properties of the PAM. 
These have apparently less influence on ferromagnetism 
in the intermediate-valence regime of the PAM. 

A quite similar picture emerges from an inspection of 
the phase diagram in figure |5| A comparison of SDA and 
MAA results shows the negative influence of quasiparticle 
damping on ferromagnetism. Furthermore, another impor- 
tant point should be noticed here: in the MAA calculation, 
a hybridization strength of = 0.2 already strongly re- 
duces the magnetic region compared to the V — Q.l case. 
It is apparent that only the region with antiparallel s-f 
coupling is affected by this. The region with parallel s-f 
coupling remains almost unchanged. It was shown in refer- 
ence that the MAA shows an anomalous hybridization 
dependence in this region. Whereas for other electron den- 
sities, an increasing hybridization strength V quickly sup- 
presses ferromagnetism]^, here the Curie temperature in- 
creases with increasing V after going through a minimum. 
The same behaviour is found within the SDA. However, 
the MPT does not show this behaviour at all. As dis- 
cussed in reference the local moments get quenched 
by formation of local Kondo singlets with increasing V . 
Being of low-energy nature, this effect is not covered by 
the other approximation schemes. The apparent stability 
of ferromagnetism around w 1 in these methods for 

large hybridization strengths seems therefore to be rather 
meaningless. 

A further difference between MAA and MPT is the 
upper critical /-level position Cf . The modifications in the 
MPT are due to the fact that the lower charge excitation 
joins with the Kondo resonance. The existence of the latter 
is neglected in the MAA. The MPT seems to be more 
reliable for determining this phase boundary. 

In figure ^ the Curie temperatures as function of the 
total electron density are plotted for the same model pa- 
rameters as used in figure ^. Again the conclusions are 
consistent: The quasiparticle damping which basically dis- 
criminates SDA and MAA, leads to a huge reduction of 
Tc- The inclusion of the low-energy physics, as done by the 
MPT, does not change Tc much, only above n^*"') k, 1.5, a 
suppression of is observed. It is further noteworthy that 
the change of sign of the conduction band magnetization 
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Fig. 7. Difference of the internal energy for the paramagnetic 
and ferromagnetic solution (see text). The inset shows the /- 
magnetization of the ferromagnetic solution. All calculations 
for (7 = 4, V = 0.1, T = and n*'"'' = 1.2. 



as seen in figure ^ does not lead to any particularity in 
the Tc curves. 

Up to now we have focused on the PAM in the 
intermediate-valence regime. What happens to the fer- 
romagnetic solution upon entering the Kondo regime? 
The situation is plotted in figure M. The inset shows the 
/-magnetization. All three methods (SDA, MAA and 
MPT) do have a self-consistent ferromagnetic solution 
for Cf well below the conduction band. An inspection 
of the internal energies (H), which can be calculated 
analogously to (0) and (^, of the ferromagnetic and 
the paramagnetic solution reveals, however, that for the 
SDA and MAA the ferromagnetic solutions are not the 
stable ones, the system is in fact paramagnetic. This is 
not the case for the MPT. Here the ferromagnetic solution 
remains stable. This indicates that in the Kondo regime, 
the low-energy properties become much more important 
concering ferromagnetic ordering. This was first proposed 
in reference where the origin of ferromagnetic order 
in the Kondo regime was identified as RKKY-like. The 
polarization of the conduction band is due to the forma- 
tion of Ko ndo screening clouds |^ . As was already argued 
in section ^.3| , the MAA and the SDA are not able to 
reproduce such a mechanism whereas the MPT should 
contain this at least qualitatively. 

So whereas in the intermediate-valence regime, the 
SDA, the MAA and the MPT show similar results, they 
give completely different pictures in the Kondo regime. 
This leads us to the conclusion that there have to be two 
distinct mechanisms driving the ferromagnetic ordering 
in these two different areas in parameter space. Whereas 
in the Kondo regime, a RRKY mechanism is doubtless 
the key factor, as discussed above and in reference the 
situation is clearly different in the intermediate-valence 
regime. Here, the correct reproduction of high-energy 
features, the charge excitation as ensured by the band- 
shift correction (llq) seems crucial. The inclusion of the 
particular low-energy properties of the PAM does not sig- 
nificantly change the behaviour. From our observations, 
we are led to propose a single-band mechanism similar 
to the one leading to ferromagnetism in the single-band 
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Hubbard model to be responsible for the ferromagnetic 
ordering in the intermediate-valence regime. 

This proposal is further supported by the following ob- 
servations: 

I) The critical interaction strength Uc is much larger in 
the intermediate- valence regime than in the Kondo regime 
thus pointing to a genuine strong-coupling effect (cf. ref- 
erence ^). 

II) The lower critical marking the breakdown of fer- 
romagnetism is in fact determined by a critical n^-^'^Q. 
The conduction electron density has no influence on the 
magnetic phase boundaries. This is a clear reference to a 
single-band mechanism. 

III) The polarization of the conduction band shows a re- 
markable behaviour as e.g. its change of sign. This, how- 
ever, does not affect important magnetic quantities as the 
Curie temperature. The polarization of the conduction 
electrons seems to be a consequence and not the cause 
of the ferromagnetic ordering of the / electrons. 

From these points, we arrive at the proposition that 
the ferromagnetic order in the intermediate-valence 
regime is due to some intra-band mechanism. For the 
Hubbard model ||lJ,Q|5^, the existence of a ferromag- 
netic phase in the strong-coupling regime was confirmed 
in the limit of infinite dimensions [I 



The mecha- 



nism driving this transition is simply based on a gain 
of kinetic energyjsj,^. This is supported by the strong 
dependency on the shape of the free {U = 0) density 
of states |5^,Q. It was confirmed that ferromagnetism 
is most favored in case of a non-symmetric DOS which 
has a divergence at or close to one of its edges. Going 
back to the PAM, we note that the hybridization leads 
in the intermediate-valence regime to an effective /-/- 
hopping. The /-electrons form a strongly correlated band. 
This band fits well into the prerequisits of a single-band 
ferromagnet as lined out above: The band is narrow, 
strongly asymetric and most of its spectral weight is, 
for appropriate values of e/, located near its edge (cf. 
figure thick lines). The proposed similarity between 
ferromagnetism in the intermediate-valence regime of the 
PAM and the Hubbard model manifests itself also in the 
fact that in both cases, the fulfillment of the high-energy 
expansion (|ll|) (as done by SDA, MAA and MPT) seems 
crucial for a proper description of the phenomenon p9| . 



of the high-energy expansion of the self-energy leads di- 
rectly to the spectral density approximation. The inclu- 
sion of quasiparticle damping effects without loosing the 
correct high-energy behaviour is possible via the modfied 
alloy analogy procedure. Finally, the modified perturba- 
tion theory still recovers the correct high-energy expan- 
sion, includes quasiparticle damping effects, and addition- 
ally incorporates, at least qualitatively correct, the special 
low-energy properties of the PAM. 

The resuhs of the SDA, MAA and MPT compare well 
in the intermediate-valence regime. By comparing with 
the Hubbard-I approximation, it becomes clear that the 
correct reproduction of the high-energy behaviour is cru- 
cial for a correct description of ferromagnetism in this pa- 
rameter regime. The infiuence of quasiparticle damping 
is, as expected, a reduction of the magnetic stability as 
indicated by a strongly reduced Curie temperature and 
smaller ferromagnetic area in the T = phase diagram. 
The low-energy physics seem to have only minor effects 
on the ferromagnetic properties. 

In the Kondo regime, the picture is completely differ- 
ent. Only the MPT yields a stable ferromagnetic phase. 
This, however, is in agreement with QMC calculations Q. 
The SDA and MAA fail to recover these results. The origin 
of the ferromagnetic ordering in the Kondo regime is an 
RKKY exchange as discussed in reference |^ . We showed 
in this paper why the SDA and MAA cannot reproduce 
such a mechanism. 

However, the good qualitative agreement between the 
SDA, MAA and MPT results in the intermediate-valence 
regime let us believe that here the driving mechanism to- 
wards the ferromagnetic transition must be of different 
nature. Our results gave some hints that this mechanism 
is similar to that driving the ferromagnetic ordering in 
the single-band Hubbard model with a band formed by 
the effectively delocalized / electrons. 

We acknowledge fruitful discussions with our coUegues A. Ra- 
makanth and G. G. Reddy from Kakatiya University, Warangal 
(India), in the process of adapting the MAA to the periodic 
Anderson model. This collaboration was financially supported 
by the Volkswagen foundation. One of the authors (D. M. ) fur- 
ther acknowledges support from the Friedrich-Naumann foun- 
dation. 



4 Summary 

In this paper, we have discussed ferromagnetism in the pe- 
riodic Anderson model (PAM), and possible mechanisms 
driving the magnetic ordering. 

We have reviewed a series of approximation schemes, 
from the Hubbard-I, via the spectral-density approxima- 
tion (SDA) and the modified alloy analogy (MAA) to the 
modified perturbation theory (MPT). This series repre- 
sents a subsequent improvement according to several ex- 
actly known properties of the model. The Hubbard-I ap- 
proximation is exact in the two limiting cases of vanishing 
hybridization {V — 0) and interaction {U = 0). Its system- 
atic improvement with respect to the correct reproduction 
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